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We investigate the non-radiative decay during nucleosyn- 
thesis of a massive tau neutrino with mass 0.1-1 MeV into an 
electron neutrino and a scalar or pseudoscalar particle, <j>. The 
full Boltzmann equation is used and shown to give markedly 
different results than the usual non-relativistic formalism for 
relativistic or semi-relativistic neutrino decays. Indeed, the 
region we investigate is where the formalism that has pre- 
viously been applied to solving this problem is expected to 
break down. We also compare the nucleosynthesis predic- 
tions from this scenario with results from the standard model 
and with some of the available observational determinations 
of the primordial abundances. It is found that for relativistic 
or semi-relativistic decays the helium abundance can be sig- 
nificantly lowered without changing other light element abun- 
dances. Since a problem with the standard model of Big Bang 
nucleosynthesis is that helium appears to be overproduced, a 
decay of the type we discuss can be a possible solution. 
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I. INTRODUCTION 

A number of indications seem to point to neutrinos 
having a mass. First of all there is the solar neutrino 
problem which is by far the strongest evidence we have 
for neutrino mass [jjj. Furthermore, there are also indica- 
tions of a non-zero neutrino mass from atmospheric neu- 
trino data 0] and, lastly, one group claims to have seen 
evidence for neutrino oscillations in a laboratory neu- 
trino beam also indicating a non-zero neutrino mass ||. 
Now, from laboratory experiments an upper limit to the 
tau neutrino mass can be obtained, which is presently 
m < 24 MeV 0. From cosmology, we have the well 
known limit on stable low mass (m < GeV) neutrinos |^] 

9u m„ 



Q u h 2 = 



2 93.03eV 



(1) 



using a present photon temperature of 2.736 K. h is the 
dimensionless Hubble constant and is the cosmological 
density parameter. g v = 2 for one flavor of neutrino and 
antineutrino. Since observations demand that fl^h 2 < 1 
H , we have a mass limit on any given stable neutrino []. 



Thus, any neutrino with mass in the range 100 eV - 24 
MeV is necessarily unstable. There are, however, many 
possible modes of decay for a massive neutrino. 

For example there is the predicted decay j7j v x — > 
Vje + e~ if the mass is larger than 2m e and the mixing 
angle between the two neutrinos is different from zero. 
A flavour changing neutral current can also lead to the 
decay Vi — > VjUjVj. There can also be other more ex- 
otic modes of decay, for example decay via emission of 
scalars or pseudo-scalars. This decay mode is generic for 
example in the majoron models of neutrino mass ||. 

The effect of such unstable tau neutrinos on nucleosyn- 
thesis have been investigated many times in the literature 
|)] |l5| , the most recent investigations being those of Do- 
delson, Gyuk and Turner (l3) and Kawasaki et al. jl4],[l5| . 
Dodclson, Gyuk and Turner have performed a detailed 
study of several possible decay modes in the context of 
non-relativistic decays, whereas Kawasaki et al. have 
performed a calculation using the full Boltzmann equa- 
tion for the decay mode v T — ► v^tfi |I^,[l7]]. In all cases 
it is found that it is possible to change significantly the 
primordial abundances via decay of the tau neutrino. 

In the present paper we focus on the decay 



(2) 



1 Note that this relation changes slightly if the heating of 
neutrinos from e + e~ annihilation is included B. 



where v T is assumed to be a Majorana particle and <f> is 
a light scalar or pseudoscalar particle. This differs from 
the decay v T — > i/^tfi in that it includes an electron neu- 
trino in the final state. Since v e enters directly into the 
weak interactions that interconvert neutrons and protons 
this decay can potentially alter the outcome of nucleosyn- 
thesis drastically. Indeed the non-relativistic results of 
Dodelson, Gyuk and Turner indicate that the primordial 
helium abundance, Yp, can be changed radically, either 
increasing or decreasing Yp depending on the mass and 
lifetime of the tau neutrino. 

Now, in the last few years, evidence has been gather- 
ing that the standard picture of the way light nuclei are 
formed in the early Universe may be facing a crisis [ jl8| . 
The main point is that helium is overproduced relative 
to the other light nuclei so that the standard theoreti- 
cal predictions are only marginally consistent with the 
observational results |M. Other measurements of the 
primordial helium abundance do yield somewhat higher 
values |H| , and the unknown systematical errors both in 
observations and in chemical evolution calculations may, 
however, be larger than presently assumed so that it is 
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perhaps premature to talk of a real "crisis" for Big Bang 
nucleosynthesis . 

Our approach will not be so much to discuss the spe- 
cific limits from nucleosynthesis since these are still quite 
uncertain as it will be a discussion of the differences be- 
tween our way of solving the Boltzmann equations and 
those previously used. Nevertheless, in light of the possi- 
bility that some new element is missing from the standard 
nucleosynthesis calculations we think that it is important 
to try and find methods of changing the nucleosynthesis 
predictions by including plausible new physics in the cal- 
culations. One possible way of doing this is to include a 
massive and unstable tau neutrino. 

In order to obtain good fits to the observational data 
it is, as just mentioned, necessary to lower the helium 
abundance somewhat compared to the other light nuclei. 
This can be achieved by having relatively low mass tau 
neutrinos decay while they are still relativistic or semi- 
relativistic. However, this is exactly the region where 
the non-relativistic formalism breaks down because it as- 
sumes a delta function momentum distribution of the de- 
cay products and neglects inverse decays. It is therefore 
of significant interest to investigate this decay using the 
full Boltzmann equation in order to calculate abundances 
in this parameter region. 

In the present paper we calculate the expected primor- 
dial abundances for a tau neutrino mass in the range 0.1 - 
1 MeV. In Section II we describe the necessary formalism 
needed for this calculation. In Section III we discuss our 
numerical results. Section IV contains a description of 
our nucleosynthesis calculations compared to the obser- 
vational data and finally Section V contains a summary 
and discussion. 



II. NECESSARY FORMALISM 

The fundamental way to describe the evolution of dif- 
ferent particle species in the early Universe is to use the 
Boltzmann equation 
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where the sum is over different possible collisional terms 
for the given particle, such as decay, scattering and pair- 
annihilation. In our case, we include the standard weak 
interactions of neutrinos with each other and with elec- 
trons and positrons. Furthermore we include a decay 
term. We shall assume, however, that the scalar par- 
ticles are collisionless except for the decays and inverse 
decays. That is, they have no self interactions and no 
other interactions with neutrinos. This may or may not 
be a good assumption, depending on the various coupling 
constants. It greatly simplifies the calculations, however. 
Now, the various terms in the Boltzmann equation can 
be written as follows 



df dm df 

Since there are only 2-particle interactions like 1 + 2 — » 
3 + 4, Cwcak can be written as 

C*wcak[/] = / d 3 p 2 d 3 p ;i d 3 p 4 A{f 1 , / 2 , h, U) x (5) 

S I M 1 12-34 5 4 ( P i +P2-P3- P 4 )(2vr) 4 . 

where A(/ l7 / 2 ,/ 3 ,/4) - (l-/i)(l-/ 2 )/ 3 /4-(l-/ 3 )(l- 
fi)fif2 is the phase space factor, including Pauli block- 
ing of the final states, and d 3 p = d 3 p/((2ir) 3 2E). S is 
a symmetrization factor of 1/2! for each pair of identi- 
cal particles in initial or final states pQ] , and \M\ 2 is 
the weak interaction matrix element squared and spin- 
summed. The matrix elements for the relevant processes 
have been compiled for example by Hannestad and Mad- 
sen pl[ . pi is the four- momentum of particle i. 

Since we are only looking at Majorana neutrinos the 
decay terms are quite simple. Since there is almost no 
net lepton number in the early Universe the Majorana 
neutrino is effectively an unpolarised species. However, 
this means that there can be no preferred direction in the 
rest frame of the parent particle. Therefore the decay is 
necessarily isotropic in this reference frame. In this case 
the decay terms can be written as |22|] 
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dE Vr A(J Vr , (8) 



where A(/„ T , /„., U) = f Vr {l - /„.)(! + U) - U.U0- ~ 
f Vr ), ml^ml T -2{ml + ml c ) + (ml-ml <; ) 2 /ml T . r is 
the lifetime of the heavy neutrino and g is the statistical 
weight of a given particle. We use g VT = g Ve = 2 and 
g$ = 1, corresponding to <f> = <f>. This assumption is not 
significant to the present investigation. Furthermore we 
shall assume that the masses of v e and <p are effectively 
zero during nucleosynthesis. 
The integration limits are 

E%{ Ei ) = ^^mi + 4( m! /m„) 2 ) 1/2 ± (9) 
-[E Vr (1 + 4(m 4 /m ) 2 ) 1 / 2 ± Pvr ] (10) 



and 
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where the index i = v e , <fi. 

Apart from the Boltzmann equation one needs equa- 
tions to relate the evolution of time, the cosmic expansion 
rate and the photon temperature. These quantities can 
be calculated by use of the energy conservation equation 

d(pR 3 )/dt+pd(R 3 )/dt = (11) 

and the Friedmann equation 

H 2 = 8irGp/3. (12) 

R is the cosmological scale factor, H is the Hubble pa- 
rameter and p is the total energy density of all particles 
present. 

III. NUMERICAL RESULTS 



no longer increases relative to that of a single standard 
massless neutrino species. This difference is clearly seen 
between different tau neutrino lifetimes. 

In Fig. 2 we show the spectral distribution of the elec- 
tron neutrino for a tau neutrino mass of 0.5 MeV and 
different lifetimes. To understand this plot better we can 
define a "relativity parameter" , p, for the decay 

m t T v T 

^ = 9MeW (13) 
A particle shifts from relativistic to non-relativistic at a 
temperature of roughly T ~ m/3. When the Universe is 
radiation dominated 



We have solved the Boltzmann equation for the evolu- 
tion of distribution functions together with the energy 
conservation equation , Eq. (|ll|), and the Friedmann 
equation, Eq. (|l2|). Specifically we have solved for masses 
of 0.1-1 MeV and lifetimes larger than 0.1 s. 
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FIG. 1. The energy density of the different neutrinos and 
the scalar particle in units of the energy density of a standard 
massless neutrino for a tau neutrino mass of 0.5 MeV. The 
full line is for r = 1 s, the dotted for r = 10 s, the dashed for 
t = 100 s and the dot-dashed for r = 1000 s. 

In Fig. 1 we show the evolution of energy density in 
neutrinos and the pseudoscalar particle for a tau neutrino 
mass of 0.5 MeV. The energy density evolves quite differ- 
ently in the different cases. Since the energy density in a 
non-relativistic species only decreases as R~ 3 compared 
to R~ 4 for relativistic particles the rest mass energy of 
the tau neutrino will dominate completely at late times if 
it is stable. If it decays the rest mass energy is transferred 
into relativistic energy so that the total energy density 
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FIG. 2. The electron neutrino distribution at asymptoti- 
cally low temperature (after complete decay) in units of the 
distribution of a standard massless neutrino. The tau neu- 
trino mass is 0.5 MeV. The full line is for r = 1 s, the dotted 
for t = 10 s, the dashed for r = 100 s and the dot-dashed for 
r = 1000 s. 

Therefore, if the decay is relativistic, 

r < t(T — m/3) * ^s. (15) 

Thus, if pi < 1 the decay is relativistic, whereas if pi > 1 
it is non-relativistic. For lifetimes of 1, 10, 100 and 1000 
s the relativity parameters are respectively 0.028, 0.28, 
2.78 and 27.8. For non-relativistic decays the decay neu- 
trino distribution assumes a rather narrow shape coming 
from the delta function energy distribution. For very 
short lifetimes the decay installs an equilibrium between 
u e , v T and (f> because of rapid inverse decays. This can 
lead to a significant depletion of high momentum electron 



neutrinos as also noted by Madsen |12| who treated this 
case of very short lifetimes using equilibrium thermody- 
namics. 
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IV. NUCLEOSYNTHESIS EFFECTS 

In order to estimate the effect on nucleosynthesis, we 
have employed the nucleosynthesis code of Kawano [p3| , 
modified in order to incorporate a decaying neutrino. 
This includes taking into account the changing energy 
density as well as the change in electron neutrino distri- 
bution. 

A decaying tau neutrino can affect nucleosynthesis in 
several different ways. Firstly, the cosmic energy density 
p is changed. Since the cosmic expansion rate is given 
directly in terms of this energy density via the Fried- 
mann equation, Eq. (|l2|), it is also changed. It is a well 
known fact that increasing the energy density leads to 
an earlier freeze-out of the n-p conversion and therefore 
produces more helium || , whereas decreasing the energy 
density decreases the helium fraction. This is the effect 
discussed by Kawasaki et al. ]l5|, namely that an MeV 
neutrino decaying into sterile daughter products while 
still relativistic or semirelativistic can actually decrease 
the cosmic energy density thereby decreasing the helium 
abundance. 

However, there is also another another effect stemming 
from the change in electron neutrino temperature. Since 
the electron neutrino enters directly into the n-p pro- 
cesses this can be called a "first order" effect and is 
potentially much more important than the "second or- 
der" effect of changing the energy density. This effect 
of change in the electron neutrino temperature has al- 
ready been discussed by several authors in the context of 
non-relativistic decays P,^3[. 
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FIG. 3. Helium abundance contours as a function of tau 
neutrino mass and lifetime for a baryon-to-photon ratio, rj, of 
3 x 10" 10 . The full line is Y P = 0.20, the dotted is Y P = 
0.22 and the dashed is Yp = 0.24. The value in the standard 
model for this r) is Y P = 0.2389. 

If the decay is non-relativistic, the energy of a produced 



electron neutrino is m/2. If this energy is significantly 
above the energy threshold for the two processes 

p + J7^^n + e + (16) 

and 

p + VI + e~ — > n, (17) 

the decaying tau neutrinos will act to produce more He 
]l3|,p[. The reason is that the absorption cross section 
at high energies is the same on neutrons and protons. 
Since there are many more protons present than neu- 
trons, more neutrons will be produced. In the end this 
leads to a higher helium fraction. However, if the mass 
of the decaying neutrino is below this threshold the pro- 
duced electron neutrinos will stimulate the conversion of 
neutrons into protons thereby actually decreasing the He 
abundance. This effect then competes with the rest mass 
effect which increases Yp. 

If the decay is relativistic the electron neutrinos are 
produced at roughly thermal energies. Effectively this 
amounts to increasing the electron neutrino temperature. 
This in turn leads to a decrease in helium production. If 
the decay takes place at high temperatures it is because 
beta equilibrium is kept for a longer time, whereas if the 
decay takes place at temperatures below the threshold 
for proton to neutron conversion it still leads to lower 
helium abundance because an increase in the electron 
neutrino temperature stimulates the conversion of neu- 
trons to protons over the inverse reaction. 

In Fig. 3 we show contour lines for the helium abun- 
dance as function of neutrino mass and lifetime. It is 
seen that helium can be significantly suppressed relative 
to the standard case if the lifetime is short enough and 
increased if the mass and lifetime are both high. If the 
mass and lifetime are both sufficiently high the helium 
abundance is instead increased. Notice also that even for 
small masses of the order 0.1 MeV the helium abundance 
can be changed significantly compared to the standard 
value for rather a large range of lifetimes. 

Our Fig. 3 should be compared with for example the 
results of Terasawa and Sato Q or Dodelson, Gyuk and 
Turner | fl3[ obtained using non-relativistic theory. The 
most straightforward comparison is with Figs. 2b and 3b 
in Ref . |9j . For long lifetimes the difference is quite small 
as would be expected since this is the non-relativistic 
limit. However, for short lifetimes the difference is signif- 
icant. For very short lifetimes, our calculated He abun- 
dance is larger than that of Terasawa and Sato. The 
reason is that if one uses the full Boltzmann equation in 
this case, decays and inverse decays will bring the parti- 
cle distributions into equilibrium as discussed in Sec. 3. 
Thus, if one keeps on going to shorter and shorter life- 
times nothing new happens since it is already decay in 
equilibrium. Therefore our curve for He flattens out in- 
stead of decreasing for very short lifetimes. For somewhat 
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longer lifetimes our He abundance is on the other hand 
smaller than that found by Terasawa and Sato. The rea- 
son here is that the decay produces a peak of very low 
momentum electron neutrinos and that these states are 
not upscattered because the weak interactions have al- 
ready frozen out. In the end this produces a somewhat 
colder electron neutrino distribution than would have 
been obtained using non-relativistic theory and therefore 
predicts less helium. This low momentum peak can be 
seen in Fig. 2 for the example of a 0.5 MeV r neutrino. 
For non-relativistic decays this peak disappears because 
low niomemtum states are not energetically accessible. 
In essence our predicted curve for the He abundance is 
therefore much flatter at short or intermediate lifetimes 
than what one would obtain using the non-relativistic 
formalism. For very long lifetimes our calculation fits 
fairly well with that obtained by Terasawa and Sato as 
could be expected. 

1 4 1 ' ' ' ' I 




osynthesis limits to tau neutrino mass and lifetime we 
will not go into too much detail regarding this point. For 
4 He we use the value calculated by Hata et al. |lj] of 

Y P = 0.232 ±0.002 ±0.005. (18) 

For deuterium the situation is somewhat complicated. 
From measurements in the local interstellar medium one 
can obtain a deuterium abundance of fl8|| 

D/H ~ 1.6 x 10~ 5 , (19) 

which can be viewed as a lower limit to the primor- 
dial abundance. However, some recent results from QSO 
absorption systems seem to indicate a primordial value 
much higher than this [ pi[ 

D/H ~ 1.9 - 2.5 x 10~ 4 . (20) 

Other similar observations yield much lower values, closer 
to the local one J25| . In light of the controversy of using 
deuterium results from these measurements, we use the 
locally obtainable lower limit in the present paper. From 
evolution arguments one can also obtain an upper limit 
to the primordial D± 3 He abundance of |2(| 

(D+ 3 Hc)/H < 1.1 x 10~ 4 . (21) 

Finally for the abundance of Li we use a bound of 

7 Li/H = 1.4 ± 0.3^Q 4 x 10~ 10 (22) 

obtained by Copi, Schramm and Turner [ p6| . 
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FIG. 4. The abundance of D, 3 He and 7 L i as a function L 
of tau neutrino lifetime. The curves have been calculated 
for m = 0.5 MeV and r\ = 3 x 10~ 10 . The full line shows ^ 1000 r 
(D/H)/1CT 5 the dashed shows ((D+ 3 He)/H)/lCT 5 and the 3 
dot-dashed shows ( 7 Li/H)/10- n . h 



In Fig. 4 we show the abundance of D, 3 He and 7 Li 
for a specific example of m = 0.5 MeV. We see that the 
abundances of these elements only change by relatively 
small amounts even for great variations in neutrino life- 
time. Thus, the main effect of the decay is to lower the 
helium abundance while leaving the other abundances 
more or less unchanged. 

The calculated abundances for different masses and 
lifetimes are compared with observational limits. Un- 
fortunately there is a great deal of controversy connected 
with these. However, since our main emphasis is on the 
differences between our approach to solving the Boltz- 
mann equations and the non-relativistic approximations 
previously used, and not so much on the specific nucle- 
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0.1 1.0 
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FIG. 5. Allowed region of tau neutrino mass and lifetime. 
The allowed region is between the two full lines. 

Altogether these are the observational values which the 
theoretical predictions should be able to reproduce. In 
the standard model the theoretical predictions are only 
marginally consistent with observations because helium 
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is overproduced compared to the other light nuclei. In 
our scenario this problem is resolved by having the tau 
neutrino decay during nucleosynthesis into an electron 
neutrino final state. 

In Fig. 5 we show the allowed region of lifetime versus 
mass for the tau neutrino using the above constraints. 
In all the mass interval from 0.1-1 MeV it is possible to 
obtain a fit to the observed abundances. Note however, 
that for masses in the high end of this region a fit can 
only be obtained in a very narrow lifetime interval. This 
is because of the very steep dependence of Yp on the 
lifetime in this region. For lower masses a good fit can 
be obtained in a broad region of lifetimes. 

Another important fact is that since the helium abun- 
dance is lowered without disturbing greatly the other 
abundances, the upper and lower bound on the baryon- 
to-photon ratio, 77, is now given essentially only by the 
limits coming from D, 3 He and 7 Li. This also means that 
a relatively high value for 77 can be accommodated, about 
6 x 10 -10 , coming from requiring that 7 Li should not be 
overproduced. 

In Fig. 6 we show the allowed region of 7710 as a function 
of tau neutrino lifetime for three different masses. We 
have also plotted the upper and lower limits to 7710 from 
the standard calculation. 
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FIG. 6. Allowed regions of 7710 = 10 10 x 77 for different 
tau neutrino masses and lifetimes. The regions inside the 
full lines are allowed regions. The vertical dotted lines show 
the consistency interval for 7710 in the standard model for our 
chosen observational constraints. 



V. CONCLUSION 



can change drastically compared to the standard value. 
This is in concordance with the findings of previous au- 
thors who used a non-relativistic treatment [p[JT^]. Our 
actual numerical values differ significantly from those 
previously obtained by use of non-relativistic formalism, 
but the general trend is the same, namely that low mass 
neutrinos decaying while relativistic or semi-relativistic 
lower the helium abundance. 

The decay we have studied differs completely from the 
Vt — * v^4> decay studied by Kawasaki et al. |l4|,|l5| be- 
cause the electron neutrinos directly affect the weak reac- 
tion rates that interconvert neutrons and protons. Only 
if much more reliable estimates of the primordial abun- 
dances are developed will it be possible to discern be- 
tween the two different decay modes. 

Our aim has mainly been to discuss the differences 
between using the full Boltzmann formalism and using 
the non-relativistic approximation in doing these calcula- 
tions. We have not done very detailed statistical analysis 
in order to obtain strict nucleosynthesis limits. 

However, it was shown that a good fit to the observed 
primordial abundances can be achieved for a large range 
of different masses and lifetimes. Given the possibly 
large unknown systematical errors in the observations 
and chemical evolution models it is perhaps too early to 
talk of a real crisis for Big Bang nucleosynthesis. How- 
ever, once the observational bounds become more strict 
there might very well turn out to be such a crisis. 

In light of this possible discrepancy between observed 
and predicted abundances we still feci that it is impor- 
tant to explore possible ways to change the light element 
abundances via plausible introduction of new physics. 
The tau neutrino decay into an electron neutrino final 
state is just such a possibility. 

Perhaps one should also finally note that even if the 
helium abundance turns out to have been significantly 
underestimated a tau neutrino decay of the type we have 
discussed can still make nucleosynthesis predictions fit 
the observations, but for completely different values of 
mass and lifetime. 

ACKNOWLEDGMENTS 

This paper has been supported by a grant from the 
Theoretical Astrophysics Center under the Danish Na- 
tional Research Foundation. I wish to thank the Max 
Planck Institute for Physics in Munich for their hospi- 
tality during the writing of this paper as well as Georg 
Raffclt and Jes Madsen who have provided constructive 
discussion and criticism. 



We have studied the decay of a relatively low mass tau 
neutrino into an electron neutrino and a scalar or pseu- 
doscalar particle using the full Boltzmann equation. It 
was found that the primordial helium abundance, Yp, 
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